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ABSTRACT 

We present the discovery of five new dwarf galaxies, Andromeda XXIII-XXVII, located in the outer 
halo of M31. These galaxies were discovered during the second year of data from the Pan- Andromeda 
Archaeological Survey (PAndAS), a photometric survey of the M31/M33 subgroup conducted with 
the MegaPrime/MegaCam wide-field camera on the Canada-France-Hawaii Telescope. The current 
PAndAS survey now provides an almost complete panoramic view of the M31 halo out to an average 
projected radius of '~ 150 kpc. Here we present for the first time the metal-poor stellar density map 
for this whole region, not only as an illustration of the discovery space for satellite galaxies, but also 
as a birds-eye view of the ongoing assembly process of an disk galaxy. 

Four of the newly discovered satellites appear as well-defined spatial over-densities of stars lying on 
the expected locus of metal-poor (—2.5 < [Fe/H] < —1.3) red giant branch stars at the distance of 
M31. The fifth over-density, And XXVH, is embedded in an extensive stream of such stars and is 
possibly the remnant of a strong tidal disruption event. Based on distance estimates from horizontal 
branch magnitudes, all five have metallicities typical of dwarf spheroidal galaxies ranging from [Fe/H] 
= -1.7 ± 0.2 to [Fe/H] = -1.9 ± 0.2 and absolute magnitudes ranging from My = -7.1 ± 0.5 to 
My = — 10.2±0.5. These five additional satellites bring the number of dwarf spheroidal galaxies in 
this region to twenty five and continue the trend whereby the brighter dwarf spheroidal satellites of 
M31 generally have much larger half-light radii than their Milky Way counterparts. 

With an extended sample of M31 satellite galaxies we also revisit the spatial distribution of this 
population and in particular we find that, within the current projected limits of the PAndAS survey 
the surface density of satellites is essentially constant out to 150 kpc. This corresponds to a radial 
density distribution of satellites varying as a result seemingly in conflict with the predictions of 
cosmological simulations. 

Subject headings: Local Group — galaxies: dwarf — galaxies: Local Group — galaxies: structure 



1. INTRODUCTION 

The last decade has seen tremendous 
understanding of our nearest giant neig 
particular highlight has been the first 
ration of the outer halo of the galaxy, 
ton Telescope wide-field imaging survey 
a ~ 100 X 100 kpc^ region down to s 
below the red giant branch (RGB) tip 



advances in our 
;hbour, M31. A 
detailed explo- 
The Isaac New- 
of M31 mapped 
» 3 magnitudes 
and revealed a 



strikingly complex pi cture of an inner halo replete with 
stellar substruct ure (llbata et al.l 120011 : [Ferguson et al.l 
I2002t llrwin et al.ll2005D. and ledt o the discovery of An- 
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dromeda XVII (|Irwin_e£iI] |2008D . With no sign of the 
substructure abating at the survey edge, attention natu- 
rally turned to the more remote regions of the M31 halo. 
Using Me gaCam on th e Canada- France-Hawaii Telescope 
(CFHT1, llbata et al.ll2007l mapped the stellar distribu- 
tion in the entire south-eastern quadrant in the radial 
range ~ 30 — 150 kpc, with an extension towards M33. 
A rich array of tidal streams and overdensities, glob- 
ular clusters and satellite galaxies were revealed. En- 
couraged by the success of these surveys we embarked 
on t he Pan- Andromeda A rchaeological Survey (PAndAS, 
McC onnachie et al.l [20091 ) . a Large Program using the 
CFHT MegaCam imager to map the entire stellar halos 
of both M31 and M33 out to distances of ^IbQ and ^50 
kpc respectively. 

As the data came in, quadrant by quadrant, in- 
creasing numbers of new dwarf spheroidal (dSph) satel- 
lite galaxies were revealed in the outskirts of M31: 
And XI, XII, XIIL XV, XVI, XVIII, XIX XX 
XXI and X XII ("Martin et al. 2006'; 'Ibata et_alj l2007t 
iMcConnachie et al. 2008; Martin et al. 200^)- In par- 
allel efforts, the existence of And IX, X and XIV was 
discovered by other groups (|Zucker et al.l 12004 120071 ; 
iMaiewski et al.l I2007D . The recently acquired second 
year of data for the PAndAS survey now yields an al- 
most complete perspective of the outer halo of M31 out 
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to '-^ISO kpc, showing for the first time the full struc- 
tural complexity of the ongoing assembly process of an 
L* disk galaxy (Fig. [Ij. Also clearly visible within 
this panoramic view are multiple dwarf galactic satel- 
lites in the M31 system including five new dwarf galax- 
ies, And XXIII, And XXIV, And XXV, And XXVI and 
And XXVII, revealed here for the first time (circled in 
red in Fig. [T]). 

However, our knowledge of the satellite system of M31 
is still limited. Current wide-area ground-based photo- 
metric observations of satellites at the distance of M31 
cannot easily reach much deeper than the horizontal 
branch, limiting detection to objects at the bright end 
of the satellite luminosity function {My ^ —6.5). If M31 
had a population of ultra-faint dwarfs like the MW, they 
would not be picked up in present surveys. Spectroscopic 
observations are currently even more constrained and are 
limited in practice to the handful of brighter red giant 
branch RGB stars close to the tip of the RGB (TRGB) 
or, if present, any bright asymptotic giant branch (AGB) 
stars. Despite these limitations much can already be said 
of the generic properties of these new objects from the 
current photometric survey data alone. 

The absolute magnitudes of the recently discovered 
M31 satel htes range from My = —6.4 for And XII and 
And XX (jMartin et al.l 120061: iMcConnachie et al.ll2008D 
to a surprisingly bright My ^ —9.7 for And XVIII 
(jMcCo nnachie et al. 200l), patently showing the incom- 
pleteness of the M31 satellite luminosity function in re- 
gions that have so far only been surveyed with photo- 
graphic plates. We note here that, even with hindsight, 
only one of these new discoveries (AndXVIII) is visible 
on the earlier photographic sky surveys and none of the 
five new satellites presented here are visible on earlier sky 
surveys either. With their relative faintness and signifi- 
cant sizes (half-light radii > lOOpc), these new systems 
are usually assumed to be (dSph) galaxies, i.e., devoid 
of any significant amount of gas and without significant 
recent star formation. The absence of gas is currently 
consistent with the results of Hi surveys but the upper 
limits on their Hi c ontent generally remains r elatively 
high, 2 - 3 X 10^ Mq (jGrcevich fc Putmanl[200l . though 
in s ome cases the H i limit is much tighter. For exam- 
ple, iChapman et al.l ()2007[) constrained the Hi mass of 
And XII to be < 3 X 10'^ Mq. In general, however, it is 
currently not possible to rule out entirely that some of 
these galaxies may still contain non-negligible amounts 
of gas. 

The properties of the surviving satellite populations 
of L, galaxies, such as M31 and the MW, are presum- 
ably directly linked to the overall evolution and assem- 
bly of the host systems. Since M31 and the MW evolved 
in a common environment, have similar morphological 
types (Sbl-II and SbcI-I I), masses (1 - 2 x lO^^ ^^^^ 
lEvans fc WilkinsonI 120001) and luminosities, one would 
expect their satellite systems also to be similar. Intrigu- 
ingly, a growi ng body of evidence suggests that this is 
not the case. IMcConnachie &: IrwinI (|2006aD compared 
the structural properties of the then known MW and 
M31 dSph populations and found systematic differences 
whereby the dSphs of M31 had half-light radii typi- 
cally twice as large as the MW dSph. A s more dSphs 
have been discovered around both galaxies. IKalirai et al.l 



(|2010D and lCollins et al.l ()2010l ) have used the kinematics 
of several M31 dSphs to show that the brighter dSphs are 
systematically dynamically colder and have lower central 
densities than their MW counterpart; while at the fainter 
end of the dSph luminosity f unctions t he structural prop- 
erties appear to overlap (Kalir ai et al .' 2010). This poten- 
tial luminosity-dependent discrepancy between the ob- 
served properties of the MW and M31 satellite systems 
is puzzling and as more satellites are discovered is worth 
reassessing. 

In this paper we analyze the properties of these new 
dwarfs and place them in the context of the rest of the 
Andromedan satellite galaxy population. Section 2 sum- 
marizes the PAndAS data and the processing steps in- 
volved. Sections 3 and 4 present the new systems and 
their structural properties and in section 5 we discuss 
them in the context of the ensemble of Andromeda satel- 
lites and compare this to the MW satellite system. 

Throughout this paper, the distance moduli of M31 is 
assumed to be (w - M)g = 24.47 ± 0.07, or 785 ± 25 
(IMcConnachie et al.ll200l . 

2. THE PANDAS SURVEY 

PAndAS builds upon the previ ous CFHT/M egaCam 
surveys of M31 pres ented in IMartin et al] (12006') . 
Ilbata et al.l (|2007D and IMcConnachie et al.l ()2(308l ) and 
we refer the reader to these papers for a full de- 
scription of the observing strategy, data reduction and 
data quality. The imaging was carried out with the 
MegaPrime/MegaCam camera mounted on the Canada- 
France-Hawau Telescope which has 36 2048 x 4612 CCDs 
with a pixel size of 0.187arcsec and 1x1 deg^ field of 
view. Although the survey is comprised of contiguous 
imaging, small gaps in exposures lead to a scientifically 
usable field-of-view of 0.96 x 0.94 deg^ for each of the 
(currently 357) pointings of the survey. Each field has 
been observed in the MegaCam g and i filters for at least 
1350 s, split into 3 x 450 s dithered sub-exposures. Good 
seeing (< 0.8") and rigorous CFHT on-site quality con- 
trol ensures that the photometry reaches g ^ 25.5 and 
i ^ 24.5 with a signal-to-noise ratio of 10. With data 
of this quality the star/galaxy separation only begins to 
degrade significantly at magnitudes fainter than this. 

Data were preprocessed by the Elixir system at 
CFHT, this includes de-biasing, flat-fielding and fringe- 
correcting the data as well as determining the pho- 
tometric zero-points. The data were then transferred 
to Cambridge where they were further processed us- 
ing a version of the Cambridge Astronomical Survey 
Unit (CASU) photometry pipeline (jlrwin fc Lewisll200l 
which was specially tailored to CFHT/MegaCam obser- 
vations. Here the astrometry of individual frames are 
refined and this information was used to register and 
then stack component images to create the final image 
products from which the survey catalogues are gener- 
ated. Finally, the astrometry is further refined, and ob- 
jects from the catalogues are morphologically classified 
(stellar, non-stellar, noise-like) before creating the final 
band-merged g and i products. The catalogues provide 
additional quality control information and the classifica- 
tion step also computes the aperture corrections required 
to place the photometry on an absolute scale. The band- 
merged catalogues for each field are then combined to 
form an overall single entry g, i catalogue for each de- 
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tected object. In this process objects lying within 1 arc- 
sec of each other are taken to be the same and the entry 
with the higher signal-to-noise measure is retained. Ob- 
jects present only on g or i are retained throughout this 
process. 

In the following, unless otherwise stated, the magni- 
tudes are presented in their natural instrumental (AB) 
system without any reddening correction (g, i). When 
required for analysis, de-reddened magnitudes {go and 
io) or reddened model isochrones have been determined 
from the Schlcgcl ct_aL (199 8) E{B~V) extinction maps, 
using the following correction coefficients: go — g — 
3.793E{B - V) and io = i - 2m6E{B - V) listed in 
their Table 6. 

3. THE PANDAS METAL-POOR STELLAR 
DENSITY MAP 

The spectacular results of the PAndAS survey are 
shown in Fig. [T] The map details the density distribu- 
tion of all candidate metal-poor stars in the stellar halos 
of M31 (M33) out to 150 kpc (50 kpc), covering more 
than 300 deg^, or around 55,000 kpc^ at the distance 
of M31. Here stars are considered as metal-poor can- 
didates if they lie in the locus of the Color-Magnitude 
Diagram (CMD) where RGB stars with metallicities of 
f« —2.5 < [Fe/H] < —1.3 fall if at the systemic distance 
of M31 ±100 kpc, and have ?-band and g-hand magni- 
tudes satisfying 21.0 < i < 24.0 and g < 25.0. In this 
we are implicitly assuming that such stars belong to pre- 
dominantly old («10 Gyr) stellar populations. 

This is the first time that a deep (Sy ~ 
33 mag/arcsec ) high resolution contiguous map of the 
majority of the extended stellar halo of any L, galaxy 
has been observed. Since the signatures of past accre- 
tions and mergers are preserved for lon ger times in the 
outer regions of massive galaxies (e.g. iJohnston et al.l 
[T996I) . the PAndAS survey has effectively provided a 
unique insight into the accretion history of M31. Al- 
though the map presented in Fig. [T] was optimized to 
search for dSph-like structures at the distance of M31, 
it also clearly reveals that M31 is looped by a series of 
giant stellar streams which have been torn off from ac- 
creted satellite galaxies, as well as other more diffuse 
stellar substructures. In particular, we note the great 
stellar arc to the North- West seen here for the first time, 
the two well-defined streams of stars crossing the South- 
East Minor axis and the diffuse debris apparently em- 
anating from M33 that is visible over a large region 
indicative of a strong ti dal interaction with M31 (see 
iMcConnachie et al.l 120091) . In addition numerous outer 
halo globular clusters are also being fou nd, many co- 
incid ent with these stream-like features (jMackev et al.l 
[201(1 . 

Interestin gly, although part of the North- West arc was 
also seen bv lTanaka et al.l ()2010[ ). their stream "F" , there 
is no evidence in our metal-poor map for any over-density 
at the location of their stream "E", nor in any of our 
other maps made with different metallicity cuts. 

Our panoramic view also offers the most complete pic- 
ture of the dwarf satellite galaxy population of M31 to 
date. Previously known dwarf spheroidals are circled in 
blue. The five new dwarf spheroidals, which are the main 
subject of this paper, are circled in red. Of the 25 cur- 



rently known Andromedan dSphs, And VI and VII lie 
well outside the field-of-view shown here, while And XIV 
lies just to the South of the surveyed region and is shown 
with a blue circle. Larger satellites of M31 such as M32 
and NGC205 are hidden within the dense inner halo, but 
NGC147 and NGC185 make an appearance at the North- 
ern edge of the image. Including M31 and M33, but not 
counting the numerous streams and other extended sub- 
structure, the number of galaxies in the Andromedan 
system comfortably visible in the surveyed area in this 
image totals 28. Since these objects are distributed all 
over the area surveyed with no obvious curtailment as a 
function of radius, it is clear that even this large-scale 
panoramic view is still not large enough to capture all 
the detail of the M31 outer halo systems. 

4. ANDROMEDAS XXIII - XXVII 

The five new dwarf galaxies were straightforward to 
detect as over-densities in the matched-filter surface den- 
sity maps of metal-poor RGB stars. There are no other 
obvious candidate satellites visible in this map or in our 
equivalent average metallicity and metal-rich maps. For 
simplicity and following convention, we have named them 
Andromeda XXIII - XXVII (And XXIII - XXVII) af- 
ter the constellation they are found in. As noted earlier 
the locations of all previously known dSphs (blue circles) 
and the five new dwarfs (red circles) are readily visible 
in Fig. [TJ Moving from East to West (left to right in 
the Figure) the new galaxies are And XXIII, And XXIV, 
And XXVII, And XXV and And XXVI. Their central 
co-ordinates are listed in Table 1. A zoom-in view of 
the distribution of stars seen in each galaxy is shown 
in the left-hand panels of Fig. [2] where the distribution 
of the candidate metal-poor RGB stars are highlighted 
with bolder dots. The other points plotted show the rest 
of the catalogued stellar population in these fields and 
serve to highlight where gaps in survey coverage might 
have an impact. The overlaid ellipses are set at twice the 
derived half-light radii and also use the ellipticity and po- 
sition angle from the likelihood analysis (see section ITT]) . 
And XXIII - XXVII all appear as fairly obvious over- 
dense concentrations of metal-poor stars in the localized 
spatial distribution maps. 

CMDs of stars lying within two half-light radii of the 
center of each galaxy (within the ellipses in Fig. [SJ listed 
in Table 1) are plotted on the left-hand side of Fig. [3] 
and are compared to CMDs of nearby reference regions 
of the same area (middle panel of Fig. [3]). These partic- 
ular reference regions were generally defined using an el- 
liptical annulus just beyond four half-light radii from the 
center of each dSph. However, for And XXVII this re- 
gion intersected too much of the enveloping North- West 
arc (stream) and an ellipse 5 arcmin to the south and 
west was used instead, though even here there is some 
contamination from stream stars. The dwarf galaxies ex- 
hibit clear RGB sequences with 0.8 < g — i < 1.5 and 
i > 21.0 which are either absent, or nothing like as ap- 
parent, in the reference regions. The narrow color width 
of the RGBs, typical of dSph galaxies, suggests that their 
stellar populations do not contain a wide range of metal- 
licities. There is no evidence that any of these galax- 
ies contains a young main-sequence population. Refer- 
ence regions are mainly devoid of RGBs and all contain a 
'red cloud' of foreground dwarf stars from the Milky Way 
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Fig. 1. — Surface density map of stars with colors and magnitudes consistent with belonging to metal-poor red giant branch populations 
at the distance of M31. The almost uniform underlying background is mainly contributed by foreground stars in the Milky Way together 
with a small residual contamination from unresolved compact background galaxies. All of the previously known M31 dwarf spheroidals in 
this region covered by the survey are readily visible as well-defined over-densities and are marked with blue circles. The five new dwarf 
spheroidals are highlighted in red. (And XIV is the dwarf spheroidal just south of the present survey area, while AndVI and AndVII lie 
respectively well to the West and North of the region shown.) NGC147 and NGC185 appear at the top of the map and M33 at the bottom 
left. The green circle lies at a projected radius of 150 kpc from the center of M31 within which most of the survey lies. In addition to the 
satellite galaxies numerous stellar streams and substructures are visible. Although the majority of small over-densities are satellite galaxies 
of M31, a few to the southern end of the map (not circled) are resolved globular cluster systems picked out surrounding nearby low redshift 
background elliptical galaxies. 



(MW) disk at the top right of the CMD with g~i> 2.0 
and i < 23.0. The fainter group of objects centered at 
0.0 < g — i < 1.5 and i > 24.0 and most readily visible in 
some of the comparison regions are generally caused by 
contamination from mis-classified compact background 
galaxies. 

The right hand panel of Fig. [3] displays the i-band lu- 
minosity functions (LFs) of And XXIII - XXVII (black) 
computed within two half-light radii of the centers, and 
a scaled LF from a nearby reference field of nine times 
larger area (red). These large area reference fields have 



also been used to measure the distances, magnitudes 
and mctallicities referred to later in this paper. For 
And XXIII, And XXIV, And XXV and And XXVI the 
reference regions are described by elliptical annuli with 
the same ellipticity and position angle of the dSph in 
question (see Table 1), the inner boundary lies at four 
half-light radii from the center of the dSph and the outer 
boundary is positioned so that the area covered is nine 
times larger (after allowing for gaps and edges) than the 
area used for the dSph. For And XXVII a circular refer- 
ence region was chosen one degree to the North-East to 
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avoid the extended tidal stream as much as possible. 

In Fig. [3] only candidate RGB stars, which were se- 
lected from the regions marked with dashed red lines on 
the dwarf galaxy CMDs, are included in the LFs. The 
z-band LFs of the corresponding reference regions, nor- 
malized from a much larger area to reduce shot noise, 
are unable to match the shape and number over-densities 
due to the dSph galaxies. At bright magnitudes, the LFs 
peter out and eventually stop beyond i ^ 21.0. This 
marks the point (TRGB) where H-shell burning ceases 
and the He-flash is triggered. The Lband luminosity of 
the TRGB is almost ind ependent of metallicit y for dSph 
systems such as these (|Bellazzini et al.l I2001D and pro- 
vides an excellent distance indicator for well-populated 
RGBs. However, since most of the new dSphs are rela- 
tively sparsely populated and the TRGB is not always 
well-defined, we have based our primary distance esti- 
mate on the luminosity of the horizontal branch and indi- 
cated the predicted location of the TRGB with an arrow 
in Fig. [3] (see Sectiong^). 

And XXHI is the largest of the newly discovered 
dwarf galaxies. Given its size and stellar density 
it may seem surprising that it has not been dis- 
covered before but even close inspection of the on- 
line digitized photographic plate sky surveys (e.g. 
http://archive.stsci.edu/cgi-bin/dss_forml shows noth- 
ing visible at this location. However, with central surface 
brightnesses of Ey^o w 27.5 (mag/arcsec^) for all of these 
new discoveries (see Section l473l and Table 1) their invis- 
ibility on previous sky surveys is not unexpected. 

The imaging of And XXV provides a good example of 
Murphy's Law and also illustrates the pitfalls of surveys 
like this with emphasis on large area coverage and survey 
speed at the expense of filling all the gaps between detec- 
tors, and also sometimes field edges (see Fig. [2]). Much of 
And XXV falls in the ^2 arcmin gap between the unbut- 
ted edges of two of the rows of 9 CCDs making up part 
of the Megacam array. The small amount of dithering in 
the (usual) 3 exposure dither sequence fills the tiny gaps 
between the long side of the detectors but does little for 
the much larger gap either side of the outer short edges 
of the middle set of the 2x9 3-edge-buttable detectors. 
It is hardly surprising that gaps and edges make an ap- 
pearance, even in this small sample of vignettes, since 
half of the surveyed region lies within w5-6 arcmin of 
either a field boundary or a large gap between rows of 
detectors. Fortunately, from a discovery viewpoint, most 
M31-distance dSphs are large enough to be found even 
in they land in the middle of a gap; while from a census 
standpoint, the impact of gaps on completeness is one 
of several factors that can readily be quantified, but is 
outside the scope of this paper. 

And XXIV and And XXVI have the most sparsely pop- 
ulated RGBs of the dSphs resulting in much noisier LFs. 
However, in both cases a metal-poor RGB sequence can 
still be clearly discerned, which is manifestly absent in 
the companion reference fields. And XXVI suffers some- 
what for lying in a field with a much shallower limiting i- 
band and g-band depth (by wO.S^^and 0.25™respectively) 
due to poorer than average seeing) compared to most of 
the other fields. 

And XXVII is a unique case among these galaxies as 
it may be directly connected to an obvious stellar stream 
which it seems to form part of. In Fig. [1] it can be 



seen to align with part of the great North- Western arc of 
stars, at least in projection. On the CMD, its stars have 
the same distribution of colors and magnitudes as stars 
selected from the arc implying that they have similar 
stellar populations and line-of-sight distances. This sug- 
gests the possibility that And XXVII is the partially dis- 
rupted remnant of a strong tidal disruption event. Since 
And XXVII is embedded in (or superimposed upon) an 
extended diffuse stream it was difficult to find a clean 
comparison field in the near vicinity and it is likely that 
the handful of stars resembling an RGB in the references 
region (Fig. ^ also belong to the stellar stream (and 
therefore, potentially. And XXVII). 

4.1. Structural Properties 

The structural parameters of And XXVIII - XXVII 
were determined by analyzing the spatial distribution of 
their resolved metal-poor candidate stars using a vari- 
ant of the m a ximum likeli hood technique developed by 
iMartin et all (|2008l I2OO90 bu t in this case bas ed on 
the Press-Schechter formalism (jSchechter fc Pre'ssl [T976'). 
The latter approach to likelihood problems such as this 
allows gaps in coverage to be directly assigned as window 
functions as an integral part of the analysis and therefore 
removes the requirement for artificially filling the gaps as 
part of an iterative solution. 

Candidate stars used in the analysis are restricted to 
those used in the construction of the metal-poor distr ibu- 
tion map shown in Fi g[Il As in IMartin et all (poM ) we 
use a simple elliptical exponential profile to describe the 
over-density of stars in the galaxy, together with a con- 
stant background level. We differ further from their pro- 
cedure by estimating the background level over a much 
larger region than used in the likelihood analysis and 
thereby fixing this parameter. These background levels 
are factored into the uncertainty estimates of the profile 
parameters but have no effect on the location parameters. 
This leaves six remaining parameters to be estimated. 
The central co-ordinates (ao, <5o), half-light radius (rh), 
ellipticity (e = 1 — b/a where a and b are the major and 
minor axes scale lengths of the system) and position an- 
gle (0; measured East from North) are derived based on 
a simple grid search starting from visually determined 
initial estimates. The grid has step sizes of 0.05 in e, 5 
in (j) and 0.05' in r^- Providing that the starting cen- 
tral coordinates lie within the half-light core of the dSph 
their exact value has no impact on the derived structural 
parameters because they are updated during the solution 
and the likelihood surface is smoothly convex around the 
maximum value. The variation of the likelihood function 
over the grid also suffices to define the likelihood surface 
around the solution poi nt for use in param eter error esti- 
mates (see for example Mart in et al.l 120081 ). Specifically, 
the uncertainties in the profile parameters are quoted as 
la errors from direct analysis of the marginalized likeli- 
hood contours from the grid search. 

The central over-density /o, is effectively a nuisance 
parameter and is readily determined, given the other pa- 
rameters, by iteratively solving the non-linear equation 
din L/dfo — 0. We note that in the Press-Schechter 
method t his is not the same as using the integral con- 
straint in IMartin et al.l (|2008D . The central co-ordinates 
were initially chosen by visual inspection of the over- 
all survey map (which is constructed with an embed- 
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Fig. 2. — Left panels: The spatial distribution of stellar sources 
around And XXIII-XXVII. Small dots represent all stars in the 
PAndAS survey and large dots correspond to stars with colors and 
magnitudes consistent with metal-poor red giant branch popula- 
tions at the distance of M31. The central ellipse in each region de- 
notes two half-light radii for each dwarf galaxy using the structural 
parameters listed in Table [T] Right panels: Background-corrected 
radial profiles of the dwarfs measured using the average stellar 
density within series of fixed elliptical annuli using the parameters 
from Table [l] Allowance is made for incompleteness due to CCD 
gaps and the edges of the survey. The error bars account for Pois- 
son counting statistics and uncertainties in the derived background 
level. The overlaid curves are derived from the tabulated model pa- 
rameters (Table [TJ but are not direct fits to the data points. 




Fig. 3. — Left and center panels: Color-magnitude diagrams 
within two half-light radii of And XXIII-XXVII and of suitable 
nearby comparison regions covering the same area after correct- 
ing for gaps in the survey coverage if required. In all cases a 
well-defined RGB is clearly visible as an over-density of stars with 
0-8 $ 5 ~ * ^3 1-5 3.nd i > 21.0 that does not have an equiva- 
lent presence in the reference CMD. The rms photometric errors 
in i and g - i as a function of magnitude are displayed in the left- 
hand side of the center panel. Right panels: The j-band luminosity 
functions of the central regions of each dwarf (black) and a suitably 
scaled much larger neighboring comparison region (red). Only stars 
within the dashed polygons were included in the LFs, equivalent 
polygons were used to select stars from the reference fields. In most 
cases the luminosity functions show a clear change in counts at the 
TRGB and the arrows mark the TRGB magnitudes predicted by 
our distance estimates (see ij 14.211 . 



ded Tangent Plane World Coordinate System) and then 
updated by using contoured isopleth maps of the cen- 
tral regions of each galaxy. From the viewpoint of their 
structural properties defined by r/i, e, 0, the precise value 
of the central coordinates (ao, So) are of secondary im- 
portance and we found in practice that we could reduce 
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the dimensionaUty of the grid search by updating these 
directly from the algebraic solutions of din L/dA£^ = 0, 
din L/d/S.ri = where and A?7 are the offsets in stan- 
dard coordinates (rotated to the best fitting ellipse coor- 
dinate frame) from the current tangent point (pio, 5o) of 
the solution (see Appendix for further details). After an 
iteration of the profile parameter grid search, the central 
coordinates are updated and the process repeated until 
a satisfactory convergence is achieved. 

To construct the radial profile for display purposes, the 
average density of stars contained within fixed elliptical 
annuli and at constant position angle, was evaluated and 
the profiles were corrected for foreground and/or back- 
ground contamination using the (much larger) normal- 
ized reference region. Due to slight gradients in field 
contaminants the error in setting the background level is 
dominated by systematics as much as by Poisson noise. 
Because of this, as noted earlier, we decoupled the back- 
ground determination by fixing the value in the maxi- 
mum likelihood estimates of the other parameters and 
then directly explored the effects of plausible changes in 
this background level by re-running the estimator with a 
small range of background values. 

The right-hand panels of Fig. [2] show the background- 
corrected radial profiles derived from the metal-poor can- 
didate spatial distributions, where the error bars include 
contributions from Poisson errors and the error in deter- 
mining the local background level. The radial distribu- 
tions shown for each galaxy were derived from average 
counts in elliptical annuli using the derived ellipticities 
and position angles, with due allowance for gaps in the 
survey coverage. We note that the exponential profile 
shown overlaid for each galaxy is not a direct fit to the 
binned data points but comes directly from the maxi- 
mum likelihood analysis outlined earlier. 

The gaps in coverage vary significantly for each galaxy, 
but even And XXV and And XXVI, the worst affected, 
are only missing ~ 10% of their stellar populations and 
still allow reasonable best fit parameters to be deter- 
mined. And XXVII presents a different challenge since 
here the background model should be more complex than 
a simple constant level because of the intersecting tidal 
stream. The simple maximum-likelihood method used 
above failed to converge in this instance because of the 
complex stream-like substructure in which And XXVII 
is embedded and we employed an alternative technique 
to derive the profile parameters. To determine e and (j) 
we constructed a contoured isopleth map of the region 
around And XXVII and computed zeroth, first and sec- 
ond moments of the surface density distribution, at a 
range of simple thresholds relative to background (e.g. 
Ilrwin et al.l 12005 1). This provided a well-defined (j) but a 
relatively poorly defined e due to the complexity of the 
stream-like distribution surrounding the core. Using the 
average values of e and (/) and the weighted center-of- 
gravity from the first moments, we then computed the 
radial elliptical profile (shown in Fig. [2]) and directly fit- 
ted an exponential profile to the inner 5 arcmin zone. 

The derived parameters for all five galaxies are listed 
in Table 1. Interestingly these new satellites span the 
complete range of parameters found previously in M31 
dSphs. And XXIII has a particularly large half-light ra- 
dius of Th — 1.04 kpc. Only one other currently known 
dSph satellite companion of M31 has a larger half-light 
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Fig. 4. — Left panels: The g-band luminosity functions of the 
central regions of each dwarf (black) together with suitably scaled 
much larger neighboring comparison regions (red). Right pan- 
els: Difference luminosity functions (black) with simple Gaussian 
plus sloping background fits to the regions around the horizontal 
branches (red). An approx ima te completeness-corrected difference 
LF, as described in Section [4.2l is plotted in gray. The g-band mag- 
nitude of the peak of the fits is indicated together with an estimate 
of the rms error in the position of the peak. 

radius; AndXIX which has = 1.7 kpc. At the other 
extreme. And XXVI is the smallest of the new galaxies 
with r/i = 0.23 kpc and also among the smallest of the 
currently known M31 dSph population. 

4.2. Distances 

Most of the newly discovered dwarf galaxies are too 
poorly populated to directly determine their distances 
using traditional TRGB methods. The absence of suffi- 
cient numbers of stars near the tip of the RGB unsurpris- 
ingly leads to unreliable distance estimates (see Fig. [3]). 
Instead we decided to use the extra depth of the g-band 
for bluer stellar populations to directly measure the lu- 
minosity of the horizontal branch (HB). This feature 
accounts for the strong over-densities sloping to fainter 
magnitudes and bluer colors at the bottom of the CMDs, 
e.g. from g — i = 1.0, i = 24.5 to g — i = 0.0, i = 25.5 
on the CMD of And XXIII (Fig. E]). The i-band depth 
cuts through this region precluding direct isochrone fit- 
ting but we can still exploit the better depth of the g- 



8 



J. C. Richardson et al. 



band directly by analyzing the (7-band luminosity func- 
tion instead. 

As a benchmark, the HB magnitude of M31 is go = 

25.2, which in the typical extinction for this region trans- 
lates to an apparent magnitude oi g — 25.4 — 25.5. Since 
most of the (7-band data we have extends down to at 
least g = 26.5 mag (and to between g = 25.5 — 26.0 
with S/N = 10), we have a good chance of resolving 
the HBs of the new dSphs. However, not surprisingly, 
the comparison regions in the CMDs show that there is 
a notable component of contamination from unresolved 
background galaxies which necessitates a careful treat- 
ment of the background correction. 

The observed g-band LFs of And XXIII - XXVII are 
shown in the left-hand panel of Fig. |4] (black histograms) 
and those of the same appropriately scaled background 
reference fields (see Section^) used for the i-band LFs are 
shown in red. Since the g-band data reaches significantly 
deeper than the i-band for such stars, only a detection in 
the (7-band is required and the number of stars contribut- 
ing to these LFs is much larger than for the equivalent 
i-band LFs. The right-hand panel of Fig. U shows the g- 
band LF of each galaxy with the background-subtracted 
(black histogram) and a best fit Gaussian (including a 
term to account for the sloping background of the LF) 
to the HB over-plotted in red. The peak of the fits and 
corresponding rms errors for each galaxy are shown in 
Fig. 21 The peak in the g-band luminosity correspond- 
ing to the HB varies from g — 25.0 mag to g — 25.7 
mag. The rms error in the g-hand magnitude over this 
range is «0.1 mag and the putative HB fits are generally 
well-defined. 

To ascertain the sensitivity of the HB estimates to 
completeness effects we have generated approximate 
completeness-corrected difference LFs using a generic 
form for the completeness as a function of point-source 
S/N. Full completeness tests are outside the scope of this 
paper but we note that in comparable datasets the com- 
pleteness level is better than 90% for lOcr fluxes and typ- 
ically drops to ~ 50% by 5a. For the galaxies presented 
here, a S/N of 5(t occurs for g-band magnitudes of 26.0 - 

26.3, where the clearest impact of this an apparent steep 
fall-off in the LF of the comparison regions. The approx- 
imate completeness-corrected difference LFs are plotted 
as gray histograms in the right-hand panels of Fig. |4l 
They show that the HB is well resolved for And XXIII, 
And XXIV, And XXV and And XXVI. For these galaxies 
the number of stars reaches a clear peak (highlighted by 
the Gaussian fit) around the HB and drops down again 
before beginning to rise as the RGB progresses to fainter 
magnitudes. The best-fit Gaussians to the HB peaks of 
the completeness-corrected LFs differed from the original 
fits by < 0.02 mag, corresponding to a maximum devi- 
ation in the final distances of 9 kpc, which is negligible 
compared to other uncertainties. 

In the case of And XXVII the completeness-corrected 
LF distorts the HB peak, although this is exacerbated by 
the difficulty of finding a clean nearby comparison region. 
Consequently here we define a lower limit only for this 
galaxy. The peak of stars at g = 25.5 mag on the bottom 
right hand panel of Fig. U marks the brightest possible 
magnitude of the HB and we use this to determine the 
lower limit on the distance. 

To calculate the distance modulus we must first esti- 



mate a standard absolute magnitude, Mg for the HB pop- 
ulations typical of these dSphs. Empirical mod els which 
estim ate the closely related My of the HB (e.g. iGrattonI 
Il998f l have been predominantly determined from glob- 
ular clusters and are shown to have a weak metallic- 
ity dependence (gradient). At the range of metallici- 
ties expected here, —2.5 < [Fe/H] < —1.5, the variation 
is less than ±0.1™which suggests a simple calibration 
based on similar dSphs which have independently deter- 
mined reliable distances from TRGB measures and/or 
HB isochrone fits from HST data. To determine the 
value of the constant we additionally measured the HB g- 
band luminosity of four relatively bright and well-studied 
dSphs (And I, And II, And HI and And XVI) which 
are covered by the PAndAS survey. The distance mod- 
uli of these well-populated dSphs, which have been de- 
termined using t he TRGB method in pre vious papers, 
are 24.43 ±0.07 (IMcConnachie et al .| [2Q0l. 24.05 ±0.06 



(iMcConnachie et al.ll2004D. 24.38 ±0 . 06 (iDa Costa et all 



|2002[ ) and 23.60 ± 0.20 (|Ibata et al.l l2007D. Since these 
galaxies contain hundreds to thousands of RGB stars the 
rms errors in these measurements are small. Another 
reason for choosing these particular galaxies is because 
they tend to lie on the MW side of M31 and hence have 
clearly-defined HBs in their g-band LFs. The formal 
rms error in locating their HB magnitudes is negligible 
(~ 0.02 — 0.03 mag) compared to the systematic errors 
due to the precise mix of their stellar populations 

By comparing to the known distance moduli of And I, 
And II, And III and And XV we find values for the 
absolute magnitude Mg of 0.89, 0.89, 0.79 and 0.67 
respectively and have therefore adopted a constant of 
Mg = 0.8 ± 0.1 in the CFHT MegaCam AB magni- 
tude system. These values are consistent, within the er- 
rors, w ith those derived for red HB stars by Che n et ahl 
(poogf ) using SDSS data for eight globular clusters with 
similar metallicities. The distances estimated by this 
method are dominated by the uncertainty of Mg , which 
in turn refiects uncertainties in the mix of blue-HB, red- 
HB, and red clump (RC) populations in these galax- 
ies. Metallicity variation is most likely a secondary as- 
pect given that all of the new objects have similar over- 
aU metallicities, [Fe/H] = -1.7 to -1.9, (see section lOl) 
to the calibrating obj ects ([Fe/H] = -1.5, -1.5, -1.9 
and -1.7 from [M cConn achie et al.ll2004l: IDa Costa et al.l 
1200a llbata et all 120071 respectively). Age effects wih 
also influence the properties of the HB and RC e.g., 
IGirardi Sz Salarii (|200i), however, given their appear- 
ance on the CMDs (Fig [3]) it is highly unlikely that the 
newly discovered dSphs contain anything other than old, 
metal-poor stellar populations. To further refine these 
distance estimates deeper targeted imaging of sufficient 
depth to allow isochrone fitting to beyond the HB mag- 
nitude is required. 

The distance moduli of the five new galaxies, after due 
allowance for line-of-sight reddening, are listed in Ta- 
ble 1. The errors in the distance moduli were calculated 
by propagating the error in Mg , the error in the measured 
location of the HB (see Fig 01 and the ±3% uncertainty 
in the calibration of the photometry. And XXIV lies 
the closest to the MW at 600 ± 33 kpc, And XXIII and 
And XXVI lie at roughly the same line-of-sight distance 
as M31, 767 ± 44 kpc and 762 ± 42 kpc respectively, 
while And XXV lies furthest from us, beyond M31, at 
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812 ± 46 kpc. The difficuhy in measuring the HB mag- 
nitude of And XXVII means we can only assign a lower 
limit to its distance of > 757 ± 45 kpc. 

As an independent consistency check we have calcu- 
lated the predicted magnitude of the TRGB in the i- 
band using our HB-based distance moduli and indicated 
them with an arrow in Fig. [3] The estimated locations 
of the TRGBs of And XXIII - XXV agree well with the 
drop off in star counts at the bright end of their LFs. 
For And XXVI and And XXVII the case is less clear 
due to the small numbers of stars, increasing the con- 
tamination from foreground stars from the overlapping 
'red cloud' region. To reliably disentangle dSph members 
from foreground stars in this part of the CMP requires 
spectroscopic followup ( e.g.. iLetarte et al.]|2OO90 . 

Further supporting evidence for the distances we have 
calculated is provided in the right-hand panels of Fig [5l 
Here theoretical isochrones have been shifted to the de- 
rived distances and compared to the CMDs of the galax- 
ies. In each case the isochrones agree favourably with 
the distribution of the RGB stars. 

4.3. Magnitudes 

The absolute magnitudes (My) were computed by 
comparing the LFs of And XXIII - XXVII to the LFs of 
And I, And II and And III. These latter dSphs are also 
part of the PAndAS survey and cover a similar range 
of metallicities and (presumably) mix of stellar popula- 
tions. Absolute magnitudes of My (Andl) = —11.8 ± 
0.1, My (Andll) = -12.6 d = 0.2 and My (Andlll) ^ 
-10. 2 ± 0.3 were taken from iMcConnachie etaTI ()2004f ) 
and IMcConnachie fc IrwinI ( 2006aD . For each dSph 
background-subtracted g-band luminosity functions were 
constructed using stars within loci such as those marked 
on the CMDs in Fig El AU LFs are corrected for CCD- 
gaps and include only stars within two half-light radii 
of the center of each galaxy using the structural pa- 
rameters listed in Ta b le 1 fo r the new dSphs and in 
IMcConnachie fc IrwinI ()2006af) for the reference objects. 
Limiting the selection to two half-light radii includes 
roughly 85% of the flux and lessens the sensitivity to 
bright interlopers, but still leaves a shot-noise dominated 
estimate at the bright end. 

We then select the portion of the LFs containing all 
stars above the w 90% completeness level {g < 25.5 mag) 
and below an upper flux limit (g > 22.5 mag) for all 
eight dSphs. Assuming the LFs are similar, the maxi- 
mum likelihood scale factors that yield the best match 
between the new dwarfs and the reference set (And I, 
And II and And III), defines the difference in magnitude 
between the galaxies. 

The variation in computed absolute magnitudes for 
each reference dSph gives some idea on the uncertainty 
introduced by assuming that the LFs of all of these galax- 
ies are the same, when they may have different underly- 
ing stellar populations. Similarly shifting the LFs by 
the maximum/minimum distance moduli allowed by our 
measurements in Section 14.2! gives a handle on the un- 
certainties introduced by the distances. As an example, 
the absolute magnitudes of And XXVI computed in this 
way range from -6.6 to -7.5 mag giving My = 7.1 ± 0.5. 
The final absolute magnitudes we quote are the mean 
of each set of three results and we estimate the error to 
be ±0.5™based on the average scatter in the computed 



sets of absolute magnitudes. As a consistency check we 
repeated the estimates using the i-band luminosity func- 
tions computed in a similar way but covering the magni- 
tude range 21.0 < i < 24.5 and found results consistent 
to within ±0.1 mag in all cases. 

The derived absolute magnitudes of the new dwarfs 
are listed in Table 1 and range from My — — 10.2±0.5 
for And XXIII to My = -7.1 ± 0.5 for And XXVI. 
These absolute magnitudes were then used in conjunc- 
tion with the derived structural parameters to directly 
estimate the central surface brightnesses, Sv.o, of the 
galaxies, also listed in Table 1. All of the new dSphs 
have similar central surface brightnesses ranging from 

Sy,o — 28.0 mag/arcsec^ for And XXIII to Ey_o = 

27.1 mag/arcsec^ for And XXV. 

4.4. Metallicities 

Metallicity Distribution Functions (MDFs) of 
And XXIII - XXVII are present ed in the lef t -hand 
panel of Fig. [5l We have used the iDotter et all ()2008f ) 
isochrones in the CFHT/MegaCam g, i system (pri- 
vate comm.) to construct the MDFs so that no filter 
transformations were required. Specifically, the data are 
compared to a grid of isochrones with age = 12 Gyrs, 
[a/Fe] = 0.0 and -2.5 < [Fe/H] < -0.5, spaced in 
0.1 dex steps, which have been reddened to match the 
data and shifted to the distance modulus we determined 
for each galaxy (the E(B-V) and (m - M)o values used 
are listed in Table 1). The right-hand panel shows how 
the resulting isochrones for metallicities of [Fe/H] = 
-2.5, -2.0 and -1.5 dex compare to the RGB sequences 
of And XXIII - XXVII. The -2.5 dex and -1.5 dex 
isochrones bound the RGB and match its curvature and 
brightness well in all cases, providing further support for 
our distance estimates. Stars within two half-light radii 
of the central co-ordinates of each galaxy were included 
in the MDFs if they satisfied 21.0 < i < 24.0 and 
were bounded by the -2.5 dex and -0.5 dex isochrones. 
Each star is assigned the metallicity of the nearest 
isochrone and a histogram of the resulting distribution 
of metallicities is constructed. The MDFs have been 
background subtracted using a suitably scaled (larger 
area) reference field satisfying the same conditions. 

We have determined the median metallicity of each 
galaxy and marked it with a gray dashed line in Fig. [5l 
Within the errors, all of the galaxies have the same me- 
dian metallicity ranging from [Fe/H] — —1.7 ±0.2 dex to 
[Fe/H] = — 1.9 ± 0.2 dex. We note that the errors due 
to uncertainties in distance modulus have little impact 
on the upper part of the RGB for such low metallicity 
systems. Re-deriving the metallicities of all five dSphs 
when the isochrones are shifted by ±0.1 mag in distance 
modulus results in change of median [Fe/H] by =f0.1 dex. 
Altering the fiducial value of [a/Fe] from 0.0 to 0.4 would 
systematically change their average [Fe/H] by -0.2 dex, 
while a ±3% photometric error in g - i would result in a 
shift of ±0.1 dex. 

The MDFs of And XXIII and And XXV, the two 
most populated galaxies, have very broad distributions 
(FWHM -^0.5 - 0.8 dex). We caution that this should 
not be interpreted too literally as a large metallicity 
spread. Although there is some evidence at brighter mag- 
nitudes for a significant spread in the locus of the RGB 
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Fig. 5, — Left panels: Metallicity distribution functions (MDFs) 
for And XXIII- XXVII derived assuming a fixed age of 12 Gyr an 
alpha-element enhanc ement of [g/ Fe] = 0.0 and using the Dart- 
mouth isochrone set l(Dotter et al.ll2008t l. Only stars within two 
half-light radii, with 21.0 < i < 24.0 and colors bluer than the 
[Fe/H] = -0.5 dex isochrone were considered. A correction for 
laackground contamination was made using a suitably normalized 
(larger) nearby comparison region and the zero-point is marked 
with a blue line. The median metallicity of each MDF is indicated 
by a dashed gray line. All of the dSphs have similar median metal- 
licities ranging from [Fe/H] = -1.7 to [Fe/H] = -1.9 ±0.2. For the 
two most populated galaxies, And XXIII and And XXV, we have 
also derived the MDF of stars with i < 23.0 (red histograms) as 
a check on the impact of photometric errors. Right panels : show 
a zoom in section of the RGBs overlaid with IDotter et aTl II2008I ) 
isochrones with age = 12 Gyrs and metallicities of [Fe/H] = -2.5, 
-2.0 and -1.5 which have been shifted to the distance modulus cal- 
culated for each galaxy. The position of the isochrones is in very 
good agreement with the magnitude range and shape of the galaxy 
RGBs. 
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Fig, 6, — Top panel: The variation of half-light radius (rjj) with 
absolute magnitude (Mv) for the five newly discovered galaxies 
presented in this paper (red circles), other M31 satellites (black 
circles) and MW satellites (triangles). The source of the data is 
referred to in the text. Middle panel: The variation of central 
surface brightness (Ey q) ^ function of My. Values of g have 
been taken from the literature where available, otherwise they have 
been estimated using r^ and Mv as described in § 5.1. Bottom 
panel: The variation of inferred central luminosity density po as a 
function of My. And XIX has a value much lower than the plot 
boundary as indicated by the downwards arrow. Note both the 
factor of 1000 variation in central luminosity (and implied baryon) 
density over the whole population and the factor 10 difference in 
this quantity between the brighter (Mv < —9.0) M31 dSphs and 
their MW counterparts. 



ues. The sparseness of the RGB for the fainter dSphs 
makes it difheuh to assess any spread in the locus though 
we note that And XXIV has a particularly narrow spread 
of derived metallicities around its median value (FWHJVI 
= 0.3 dex), and is perhaps the galaxy with the clearest 
sign of having a simple stellar population. The JVIDFs 
of And XXVI and And XXVII are too noisy to make 
comparable inferences from their distributions. 



for And XXIII and And XXV, for the more metal-poor 
isochrones, photometric errors, particularly at the faint 
end, contribute significantly to the spread in the derived 
JVIDF and deeper targeted observations are needed to fur- 
ther investigate this. The red histograms in Fig.[S]are the 
IVIDFs derived for And XXIII and And XXV when only 
stars brighter than i = 23.0 are considered. Both dis- 
tributions, particularly that of And XXV, are somewhat 
narrower than before, but although photometric errors at 
fainter magnitudes contribute to broadening the IVIDF, 
they have little effect on the derived median [Fe/H] val- 



5. DISCUSSION 

5.1. Structural properties 

The five new dSphs have characteristics typical of the 
range of properties of the previously known 1V131 dSph 
population. The top panel of Fig. [5] shows their loca- 
tion in the (semi-major axis) half-light radius -v- abso- 
lute magnitude plane together with the location of the 
other 1VI31 dSph satellites and equivalent MW satellites 
satisfying My < —4.5. The values for the MW satel- 
lites Forna x, Leo I, Sculptor, Leo I I and Sextans were 
taken from llrwin fc Hatzidimitriou I fl995.) while the pa- 
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TABLE 1 

Derived properties of the satellites 



Parameter 


And XXIII 


And XXIV 


And XXV 


And XXVI 


And XXVII 


a (J2000) 


01 29™21.8'^ ± 0.5'' 


01 18™30.0'' ± 0.5'' 


00 30™08.9'' ± 0.5^ 


00 23™45.6'' ± 0.5^ 


00 37™ 27.1^ ± 0.5^^ 


5 {J2000) 


38°43'08" ± 10" 


46°21'58" ± 10" 


46°51'07" ± 10" 


47°54'58" ± 10" 


45°23'13" ± 10" 


('.i') (°) 


(131.0, —23.6) 


(127.8, —16.3) 


(119.2, —15.9) 


(118.1, —14.7) 


(120.4, —17.4) 


E(B - 


0.066 


0.083 


0.101 


0.110 


0.080 


(m - M)J| 


24.43 ± 0.13 


23.89 ± 0.12 


24.55 ± 0.12 


24.41 ± 0.12 


24.59 ± 0.12 


D (kpc) 


767 ± 44 


600 ± 33 


812 ± 46 


762 ± 42 


> 757 ± 45 


J-MSi (kpc) 


~ 126 ± 44 


~ 197 ± 33 


~ 97 ± 47 


~ 101 ± 42 


~ 86 ± 48 


[Fe/H]('^' 


-1.8 ±0.2 


-1.8 ±0.2 


-1.8 ±0.2 


-1.9 ±0.2 


-1.7 ±0.2 


rji(arcmin) 


4.6 ±0.2 


2.1 ± 0.1 


3.0 ±0.2 


1.0 ±0.1 


1.8 ±0.3: 


rhlpc) 


1035 ± 50 


378 ± 20 


732 ± 60 


230 ± 20 


455 ± 80: 


4> (N to E) (°) 


138 ± 5 


5 ±10 


170 ± 10 


145 ± 10 


150 ± 10: 


e = l-b/a 


0.40 ± 0.05 


0.25 ± 0.05 


0.25 ±0.05 


0.25 ± 0.05 


0.4 ±0.2: 


Mv 


-10.2 ± 0.5 


-7.6 ±0.5 


-9.7 ±0.5 


-7.1 ±0.5 


-7.9 ±0.5 


"^Vfi (mag/arcsec-^) 


28.0 ±0.5 


27.8 ±0.5 


27.1 ±0.5 


27.4 ±0.5 


27.6 ±0.5 



''Notes. The errors quoted for the profile parameters derived in ii l4.1l {a. S.Vh.d) and e) are the la errors from the marginahzed likeUhood 
contours of the grid search. Structural parameters derived for And XXVII are marked ':' to refiect the additional uncertainties caused 
by the presence of the intersection stellar stream. ^^^E{B — V) interpolated from the Schlegel et al. (1998) dust maps. Due to the 
sparsity of the upper RGB, distance moduli are calculated from the magnitude of the horizontal branch from the deeper g-band data. ('^^ 
The median metallicity of the metallicity distribution functions from Fig. [5] 



rame ters for the res t of th e MW satelhtes were taken 
from lOkamoto et al.l ()2008[) and Okamoto (PhD thesis) 
and include the MW ultra faint dwarfs (UFDs) Boo I, 
Hercules, CVn I, CVn II, Leo IV, Leo T and Uma 1. 
The equivalent dat a for the p reviously known M31 satel- 
lites is from McConnachie & Irwin (2006a); Marti n et a,l.l 
(l2006l): ilba ta et all ()2007[ ): iMcConnachie et al.l (|2008[ ): 
iMartin et al . (2009i). 

As expected, there is a general correlation between the 
size of the galaxy and its magnitude, though we caution 
that fainter than My « — 8 observational selection effects 
may be limiting the detectability of larger half-light radii 
systems. The MW satellite population extends to fainter 
magnitudes than M31 satellites because, at the distance 
of M31, even relatively compact dwarf galaxies fainter 
than My ~ —6 are beyond our current detection limits. 
Although this rules out any direct constraints on putative 
UFDs around M31 it is interesting that the number of 
M31 satellites per absolute magnitude interval shows no 
sign of tailing off towards the faint end. This suggests 
that an equivalent population of UFDs are waiting to be 
discovered. 

Following their analysis of the stru c tural p roper- 
ties of And I- VII, IMcConnachie fc IrwinI ()2006al) noted 
that the dwarf satellites of M31 tended to have 
larger radii than MW dwarf satellites and were typ- 
ically twice as big. Subsequ ently, 18 new M3 1 
dSphs (A nd XI - And XII I (|Martin et al.l l2006l ). 
And XIV (IMaiewski et al.l [200l. An d XV and And XV I 
(jlbata et al.l 120071 ). A^ XVII (llrwin et al. ' 
And XVIII - And XX (Mc Connachie et al 



2008) 



20081) 



And XXI and And XXII (Martin et all I2009D ) and 10 
new MW dSphs, including six with My < — 4.5, have 
been discovered: CVn I (jZucker et al.l 120061 ): Boo I 
([Belokurov et al. 2006); Hercules, C Vn II and Leo IV 
Belokurov et al. 2007) ; and Leo T (jlrwin et al.l 120071 ) 



Recent work by Kalirai et al.l ()2010[ ) found that there is 
significant overlap between the M31 and MW satellite 
populations at the low luminosity end {L < 10^ Lq). 
Fig. [S] includes data from the new M31 and MW satel- 



lite galaxies more luminous than My = —4.5 in addi- 
tion to the previously known satellites and the newly 
discovered galaxies presented for the first time in this 
paper. And XXIII-XXVII. When the new data are in- 
cluded there is no clear difference in the My — rh relation 
between MW and M31 satellites for absolute magnitudes 
fainter than My w —9. However, for the subset of dSph 
galaxies brighter than My « —9 it remains the case that 
the M31 dSph satellites are generally twice as extended 
as their MW counterparts, even though they are found 
at a similar range of distances (« 100 — 300 kpc) from 
their host galaxy. Of the new dSphs, the brighter pair. 
And XXIII and And XXV, continue this trend, whilst 
the three with My > -9 (And XXIV, And XXVI and 
And XXVII) have half-light radii typical of both M31 and 
MW dSphs in this magnitude range. It is unlikely that 
this is a completeness issue since MW satellite galaxies 
brighter than My = —9 and with rh > 600 pc should be 
easy to find if they lie w ithin the expected di stance range 
of MW satellites ( e.g.. iKoposov et al.|[2008l) . 

As an alternative way of looking at these distributions, 
the variation of central surface brightness, Sv,Oj as a 
function of half-light radius is shown in the middle panel 
of Fig. H Fainter than My « -9, MW and M31 satellites 
show a similar range of central surface brightness. How- 
ever, for those brighter than My ?» —9, M31 satellites 
have significantly fainter central surface brightness than 
their MW counterparts by an average of r^2 mag/ arcsec^. 

Finally, in the third panel of Fig. [5] we make the plausi- 
ble assumption that, to first order, we can deconvolve the 
projected surface den sity distribution by approximating 
it as a Plummer law (jPlummer II1911I) with the geomet- 
ric mean half-light radius deduced from the parameters 
in Table 1 . This allows us to analytically infer the equiv- 
alent central luminosity density and also, by effectively 
assuming a constant baryonic mass:to:light ratio, gain a 
measure of the variation of the central baryon density. 

In M31 it appears that the brighter (classical) dSphs 
are not only twice as extended as their MW counterparts 
but also correspondingly have central stellar densities a 
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factor of ten smaller. This is a significant difference and 
suggests that the local environment around these two L» 
galaxies has exerted a strong infiuence on the formation 
and evolution of these dSph systems. 

In a study attempting to relate the observed prop- 
ertie s of Local Group dSphs to their dark matter con- 
tent, |Penarrubia^eF^I7] (|2008D suggested that, under the 
assumption that M31 and MW satellites have similar 
dark matter halos, the systematic difference in their sizes 
should be complemented by a systematic difference in 
their kinematics. Specifically, they predicted that M31 
satellites should have a velocity dispersion ~ 50% — 100% 
larger than the corresponding MW satellites. However, 
a recent kinematic survey of And I, And II, And III , 
And VII, And X a nd And XIV bylKalirai et all (|20lol ) 
and similar work by iCollins et al.l (j2010f ) has shown that 
this prediction is not borne out by observation. M31 
satellites with measured velocity dispersions range from 
And X at Mv = -8.1 and = 3.9±1.2k ms~^ to And I 
at M y = -11.8 and cr„ = 10.6 ± l.lkms"^ (jKalirai et al.l 
I2010D . Meanwhile, MW satellites covering the same ab- 
solute magnitude interval have = 9.5 ± 2.0kms~^ 
(Draco at My = -8.3) to ct„ = 1 0.5 ib 2. 0kms~^ (For- 
nax at My = -13.0) acc ording tolMateo I (Il99 8). This 
suggests, as also noted bv lPenarrubia et al. I (|2608V that 
the nature of the dark matter halos themselves may be 
responsible for the observed differences between the two 
satellite systems. 

5.2. Satellite spatial distribution 

An interesting peculiarity in the M3 1 satelli te distri- 
bution noted bv iMcConnachie fc IrwinI ()2006b[ ) was the 
tendency for the then known satellites of M31 to be pref- 
erentially located on the MW side of M31 (14 out of 16). 
With the discovery of 18 more M31 satellites since then 
it is worth revisiting the satellite distribution again. As 
before, we define an M31-centric coordinate system such 
that the disk of M31 defines the galactic plane i.e., lati- 
tude 5j\/3i = 0°, and such that the MW lies at longitude 
Ij^i^^ = 180° and latitude hMzi = -12.5°. Within this 
coordinate system the distribution of M31 satellites, as 
would be seen from the center of M31, is shown in the 
Aitoff projection of Fig. [71 The outlined Great Circle 
has a pole centered on the MW position and splits the 
M31 sky into those objects lying on the MW side within, 
and on the opposite side to the MW, without. Although 
there is still an apparent bias, with 24 out of 34 satellites 
lying on the MW side of M31, the offset in the line-of- 
sight distance distributi on is much less significant tha n 
that shown in figure 8. of IMcConnachie fc IrwiJ (|2006bD . 
With the additional satellites, and updated information 
for the rest, the median offset toward the MW changes 
from 40 kpc to 21 kpc, and changes from a 3cr bias to 
less than a 2a event. 

On the other hand, the overall spatial distribution of 
the satellites of M31 reveals an unexpected result. The 
projected (surface) spatial distribution of satellites de- 
picted in Fig. [8] highlights the visual impression com- 
mented on previously from the two-dimensional overview 
of satellite locations seen in Fig.[TJ Until the main limit of 
the current PAndAS survey is reached at a projected dis- 
tance of «150 kpc from M31, there is an essentially uni- 
form surface density of satellite galaxies as first pointed 



out bv IMcConnachie et all ()2009D . The turnover near 
150 kpc in the projected cumulative distribution corre- 
sponds precisely to the limit of the main PAndAS survey. 
Since we know from their three-dimensional distribution 
that M31 satellites exist at distances up to at least 300- 
400 kpc from M31, it is quite likely that even with the 
current PAndAS survey we are still seriously incomplete 
in our census of even the relatively bright (My ;$ — 7) 
satellites. 

In compiling this figure we have included all suspected 
Andromedan satellites irrespective of their location rel- 
ative to our currently surveyed region. The inclusion, 
or otherwise, of marginal Andromedan members such 
as And XVIII, makes little difference to the conclu- 
sion. The surprising result is that out to our approxi- 
mate survey limit at sal50 kpc in projection, the surface 
density of satellites is essentially constant. This corre- 
sponds to a three-dimensional radial density distribution, 
p{r) oc r~^, a result seemingly in conflict with cosmolog- 
ical simulations. As a particular example of one of these, 
we show in Fig. |8] a standard cosmological model predic- 
tion of sub-halo density profiles based on an Einasto ra- 
dial density distribution with parameters a = 0.678 and 
r-2 = 200 kpc, taken from Figure 11 of iSpringel et al.l 
(|2008f ). Since the normalization here is arbitrary, both 
overlaid models were normalized to have the same mean 
projected surface density as observed out to 100 kpc. A 
constant surface density to sal 50 kpc in projection pro- 
vides a reasonable description of the observed satellite 
population whereas the cosmologically motivated model 
has a much steeper radial density fall-off and a cor- 
respondingly much flatter predicted cumulative surface 
density distribution. 

The present observed distribution of dwarf satellite 
galaxies will undoubtedly have been influenced by their 
range of orbital properties and evolutionary histories, 
with tidal effects preferentially destroying systems closer 
to the nucleus of M31. While we might expect that the 
original distribution of recognizable satellites could have 
been more centrally concentrated, the evolution over 
time of the distribution, and properties, of the surviv- 
ing systems is less obvious. Although detailed simula- 
tions of the distribution of dark matter sub-halos have 
been made, the problem of accurately linking the sur- 
viving sub-hal os with the observed satellite systems still 
remains ( e.g.. lLudlow"eraIll20T0l ). 

6. SUMMARY 

We have presented for the flrst time a panoramic view 
of metal-poor stars in the M31 halo out to an average pro- 
jected radius of ~ 150 kpc using data obtained with the 
MegaPrime/MegaCam wide-field camera on the CFHT 
during the PAndAS survey. Dwarf satellite galaxies and 
tidal debris streams have been discovered out to the edge 
of the survey suggesting that much more is to be found 
at still larger radii. The degree of sub-structure visible 
in the overall map highlights the continuing growth and 
evolution of the outer halo and disk of M31 and yields an 
insight into the accretion process of an disk galaxy. 

In the main part of the paper we have characterized 
the five most recently discovered dwarf spheroidal galax- 
ies uncovered by PAndAS, Andromeda XXIII - XXVII. 
They all contain stellar populations typical of dSph 
galaxies, with no obvious sign of substructure, a rela- 
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Fig. 7. — An Aitoff projection of the all-sky distribution of the M31 satellite system as viewed from the center of M31. The co-ordinate 
system is defined using Zm31 measured around the disk of M31 and locating the MW at /m31 = 0° and femsi = —12.5° to account for the 
inclination of the M31 disk along the MW line-of-sight. The dot-dash Great Circle centered on the MW splits the sky into MW near-side 
and far-side. With the extra M 31 satellites now known the earlier apparently significant asymmetry in the satellite distribution noted by 
IMcConnachie & Irwiiil I l2006bl') has now largely disappeared. 




50 100 150 200 



Projected distance (R kpc) 

Fig. 8. — The cumulative distribution of the number of satellite 
galaxies as a function of projected distance from M31. The overlaid 
models are: a power law distribution of the form p{r) oc r~" 
with a = 1 (dashed line); and an Ei nasto density distribu tion with 
a = 0.678 and r_2 = 200 kpc ( e.g.. ISpringel et al.ll200g : dot-dash 
line). Note particularly that the power law model corresponds to 
a constant surface density of satellites and that the clear break in 
the projected cumulative profile at the main PAndAS survey limit 
of f5;150 kpc suggests a large number of satellites remain to be 
discovered further out. The cosmologically motivated model does 
not provide a good fit to the observed surface density distribution. 



tively narrow RGB and have mean metalhcities ranging 
from [Fe/H] = -1.7±0.2 dex to [Fe/H] = -1.9±0.2 dex. 
Although the central surface brightnesses of the galaxies 
are broadly similar {^v,o = 27.5 ± 0.5 mag/arcsec^, they 
vary in absolute magnitude (My — —7.1 ± 0.5 to -10. 2± 



0.5) and have very different scale lengths (Rh = 230 pc 
to 1035 pc). These additional discoveries continue the 
trend whereby the brighter M31 dSphs have significantly 
larger half-light radii than their MW counterparts. 

• And XXIII is a bright (My = -10.2) satehite 
galaxy lying at the same line-of-sight distance as 
its host and ^ 126 kpc away from it. With a half- 
light radius of rh = 1035 pc it is one of the largest of 
M31's presently known satellites and almost twice 
as extended as the largest known MW dSph satel- 
lite ( Fornax: rh = 636 pc, llrwin fc Hatzidimitriou I 

[T99a . 

• And XXIV and And XXVI. These two galax- 
ies resemble the traditional dSphs belonging to the 
MW group in many ways. They have small half- 
light-radii (rh = 357 pc and 232 pc) and are rela- 
tively faint (My = -7.6±0.5 and -7.1±0.5). Both 
galaxies are very sparsely populated (see Figs[5]and 
131) • Our line-of-sight distance measurements show 
that both galaxies lie on the near side of M31 rela- 
tive to the MW and, at a distance of 600 ± 33 kpc. 
And XXIV is one of the closest M31 dSphs to us. 

• And XXV is characteristic of M31's dSph satellite 
system in that it is large (rh — 693 pc) and bright 
(My — —9.7). It has the highest central surface 
brightness of the newly discovered galaxies (Sy,o = 
27.1 mag/arcsec^), still not bright enough to have 
been found by previous shallower surveys such as 
the SDSS. 

• And XXVII appears to be in the process of be- 
ing tidally disrupted by M31. It is embedded in 
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a great Northern arc of stars which loops around 
M31. The stream of stars surrounding the galaxy 
has made it very difficult to extract a clean ref- 
erence field to characterize the properties of the 
foreground and/or background. Of the five galax- 
ies presented in this paper, And XXVII is furthest 
from the MW (^827 ± 47 kpc) and the nearest to 
M31 (Rm3i ~86 kpc). 

Finally, we revisited the spatial distribution of the 
overall M31 satellite population and demonstrated that 
the apparent bias in the satellite location with res pect 
to the MW noted by iMcConnachie fc IrwinI ()2006b[ ) has 
gradually disappeared with the discovery of more satel- 
lites. However, in recompense, the newly discovered 
satellites now allow an almost complete census out to 
a projected distance of wl50 kpc and show a surprising 
uniformity of their surface density out to this limit. This 
not only suggests that a large fraction of the M31 satel- 



lite population remains to be discovered but also that 
the three-dimensional radial density distribution of M31 
satellites varies approximately as p{r) oc out to at 
least 150 kpc. It will be enlightening to compare this 
result with more detailed cosmological simulations based 
on the estimated numbers of surviving satellites of 
galaxies. Surveying to comparable depths to beyond the 
virial radius of M31 (w250 kpc) is probably impractical 
with MegaCam but well within the capab ilities of the up- 
coming Hyper-Suprime Cam on Subaru (jMivazaki et al.l 
2006) and would correspondingly provide a more exact- 
ing test of cosmological predictions. 



We thank Sidney Van Den Bergh for a careful reading 
of the manuscript and for his helpful comments on the 
content. Thanks also to the entire staff at CFHT for 
their great efforts and continuing support throughout the 
PAndAS project. 



APPENDIX 
FITTING SURFACE PROFILES 

We want to estimate the likelihood of observing the data points at position vectors {rj} , i = 1,2, m, where 

/(r) is the surface density profile of a given dwarf galaxy superimposed on a uniform local background. Following the 
Press-Schechter formalism (Schechter & Press 1976), imagine partitioning the observed region into "cells" of surface 
area 5S. Let the expected number of observed data points in cell i be (j)i where 



</'^ = /(ri) 5S 

then the probability of observing Xi points in cell i is given by a Poisson distribution 
and therefore the hkehhood function for the ensemble of data points in the region is 



(Al) 
(A2) 

(A3) 



Let 6S 0, then Xi — 1 if a. point is observed at this location, or Xi = Q if none are detected. The likelihood function 
now becomes 

L = l[e ^ . e-^' (A4) 



where the first product is over empty cells and the second is over occupied cells. Taking natural logarithms this further 
simplifies to 



In L = — 0j + In 



(A5) 



where the first summation is now over all cells and the second only occupied cells. Since 5S Q we can replace the 
first summation with a surface integral over the region such that 



In L = 



P ni 

Jj{r) dS + lnf{v,) 



(A6) 



Unobservable regions due to gaps etc., are trivially included in the normalizing integral by use of a window function, 
which of course has no impact on the discrete summation of observed data points. 

To illustrate the use of this method, we model the dSph as an ellipse with ellipticity e, position angle in a standard 
coordinate grid, and use an exponential profile of scale length a along the major axis. In this case 



/(ri) = /„ e--^/'^ + b 



(A7) 



where using standard coordinates and A77 to denote offsets from the center of the dSph i.e., ^ — and rj — rjo^ we 
have 

rl = (A^, cosO + A77, sinOf -f (-A^i sinO + Aj]i cos9f/{l - ef (A8) 
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If the area for analysis, A, covers sufficient ellipse scale lengths, and there are no gaps in the coverage, then 

m 

lnL= -2-Kfo a^(l - e) - + ^ In {fo e"'''/" + b) (A9) 

i=l 

and for a given trial grid solution for a,e,9, and for a particular background b, the central density fo is defined by 

= = -27ra2 1-e +y -. AlO 

dfo ~{ fo e-''*/" + b 

which has a straightforward iterative solution for fo- Likewise the central coordinates of the dSph, in a coordinate 
system rotated for convenience to lie along the major and minor axes, 77^, are defined by 

^ - - V ^ /° ^""'^^ (All) 

which implies that 

{-[ fo e-'-*/« + b n ' fo e-'-*/« + 6 n ^ ' 

i.e., a simple weighted sum over the coordinates, with a similar solution for 77^. 

For numerical stability in solving the latter equation, we introduce a softening parameter in the core of the dSphs 
r? — )■ r? + (? where c is fixed to be equivalent to a 20pc scale length. This also makes some allowance for the often 
observed flattening of the inner profile of dSphs with respect to a pure exponential. 
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